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We report the observation of quantized translational and rotational motion of molecular hydrogen 
inside the cages of Ceo- Narrow infrared absorption lines at the temperature of 6K correspond to 
vibrational excitations in combination with translational and rotational excitations and show well- 
resolved splittings due to the coupling between translational and rotational modes of the endohedral 
H2 molecule. A theoretical model shows that H2 inside Ceo is a three-dimensional quantum rotor 
moving in a nearly spherical potential. The theory provides both the frequencies and the intensities 
of the observed infrared transitions. Good agreement with the experimental results is obtained by 
fitting a small number of empirical parameters to describe the confining potential, as well as the 
ortho to para ratio. 
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Endohedral complexes of H2 molecules trapped inside 
fuUerene cages have been synthesized recently [l|, y, Q- 
Apart from their chemical interest and importance, these 
remarkable systems are ideal testbeds for the study of 
diatomic quantum rotors in a confined environment. 
H2@C6o is different from quantum rotors studied so far, 
which were two-dimensional or showed hindered rotation, 
like H2 on a Cu surface[j|, H2 in intercalated graphite 5] 
or C2 in a metallofullerene[6[. The small mass and large 
rotational constant of H2 makes it the least sensitive of 
molecules to the corrugations of the potential surface. 
Also, Ceo provides a nearly spherical bounding potential. 
Theoretically it has been shown that in Ceo cages the 
quantum rotors C0[7| and H2[8|,|9|] should have a measur- 
able translation-rotation coupling, a feature that has not 
been experimentally resolved to our knowledge. More- 
over, there is very little experimental information on the 
quantum dynamics of H2 in fuUerene cages Ifl lul llJ, ll3| . 

In general, isolated homonuclear diatomics have no in- 
frared (IR) activity [ij]. However, H2 does display IR 
activity in situations where there are intermolecular in- 
teractions present, such as in the solid and liquid phases 
ISl lld | , in constrained environments IT], [iSJ, |l9|, |20| , and 
in pressurized gasses[2ll |22|- IR spectra of such systems 
are usually broad due to inhomogeneities in the system 
or due to random molecular collisions. As an exception, 
narrow lines are observed in semiconductor crystals [23J 
and solid hydrogen [24J. Similarly we expect narrow lines 
in solid H2@Ceo, where the broadening of IR lines is sup- 
pressed by homogeneous distribution of trapping poten- 
tials provided by Cgo molecules and by weak van der 
Waals interactions between the molecules. 

In this Letter we study the dynamics of H2 in cages 
of Ceo in the solid state with infrared spectroscopy. The 
observed spectra are described by a three-dimensional 



quantum rotor confined in a nearly spherical potential 
exhibiting translation-rotation coupling. 

The H2@Cfio powder sample (10 mg) was prepared as 
described in 3] and pressed into a d = 0.25 mm thick 
pellet. IR transmission measurements were made with 
an interferometer Vertex 80v (Bruker), halogen lamp, 
and MCT detector with an apodized resolution 0.3 cm~^. 
The sample and the reference open hole were inside an 
optical cryostat with KBr windows. The absorption coef- 
ficient a{uj) was calculated from the transmission Tr{uj) 
through a{uj) = -d-'^\ii[TriLu){l - R)-'^] with the re- 
flection coefficient R = [{n — l)/(?i +1)] . A frequency- 
independent index of refraction n — 2 was assumed [25| . 
The low temperature IR absorption peaks of Il2@C6o are 
located in four narrow spectral bands between 4060 and 
4810 cm~^, see Fig.[T] This region corresponds to the 
H2 stretching mode and its rotational/translational side- 
bands. Peaks in 4250 cm~\ 4600 cm"^ and 4800 cm~^ 
regions (panels b, c, and d in Fig.[T|) are assigned to 
vibrational excitations of II2 accompanied by transla- 
tional and/or rotational excitations. It is the translation- 
rotation coupling that splits the (5(1) line, shown in 
Fig-Uti into three peaks. Weak transitions around 
4070 cm~^ (Fig-Uti-) represent pure vibrational excita- 
tions of the II2 molecule and are forbidden in the ap- 
proximate theory presented below. 

The position and orientation of the II2 molecule is de- 
scribed using spherical coordinates R = {R, Q, $} and 
r = {r, 6, (j)} where R is the vector from the center of the 
Ceo cage to the center of mass of H2 and r is the inter- 
nuclear H-H vector. As first approximation we consider 
decoupled translational, rotational and vibrational move- 
ment of II2. The translation of the confined molecule 
may be treated using the isotropic three-dimensional har- 
monic oscillator model. The appropriate translational 
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FIG. 1: Baseline-corrected IR absorption spectra of H2@C6o 
at 6 K (black) and the best fit theoretical spectrum (grey) in 
panels (b), (c), and (d). All transitions are from the vibra- 
tional state u = to i; = f . The letter indicates the change 
in J (Q for AJ = and S for AJ = 2) and the number in 
parentheses is the initial J value. For the assignment inside 
a rotational branch, see Table |T] and Fig. [5] (a) Fundamental 
vibrational transitions; these are forbidden in the first-order 
theory so a simulated spectrum is not shown, (b) AJ — 0, 
A7V — +1 transitions, (c) para-H2 and (d) ori/io-H2 S tran- 
sitions with AN = -1-1, AJ = -1-2. 



quantum numbers are N = 0, 1, . . ., the orbital angular 
momentum quantum number L, which is an integer with 
the same parity as N and the azimuthal quantum num- 
ber Ml. The radial part of the wavefunction depends 
both on N and L. The translational eigenfunctions are 
*I^lmJ^:0.^) = *?;i(i?)lLA/je,$) where the radial 
wave function yXfr and the spherical harmonics Ylml 
are defined in [26|. The rotational wavefunctions, de- 
fined by the rotational quantum numbers J — 0, 1, . . . 
and Mj = — J, — J+l,...,+J, are given by the spheri- 
cal harmonics Yjmj (^, 0)- It is convenient to use bipolar 
spherical harmonics with overall spherical rank A and 
component Ma, defined as follows: 



FkLm 



E Ct^ltjMYLMAQ,'^)YjM,(e,<l>), (1) 
Ml.Mj 

- |9, $, 6', 0} and C are the Clebsch-Gordan co- 
27| . The full wavefunction describing the mo- 



where Vi - 
efficients 

tion of the H2 molecule may be written as \vJNLAMa) — 
*^(r)4'^i(i?)FfJ^^^(0) where *^(r) is the vibrational 
wavefunction with a quantum number v. The total nu- 
clear spin / of the H2 molecule determines whether it is 
either in a para state (7 = and J even) or in an or- 
tho state (7 = 1 and J odd). The Hamiltonian H for the 
trapped molecule includes coupling terms between the vi- 



brational, translational, and rotational motion. For sim- 
plicity, we neglect all matrix elements non-diagonal in v 
and introduce a parametric dependence on v. 
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where ^Ji^R is the vibration-rotation Hamiltonian, p is 
the molecular momentum operator and m is the molecu- 
lar mass. The superscript prefix v is used to indicate an 
implied dependence on the vibrational quantum number. 
■"T^ is the potential energy of a molecule at a given posi- 
tion and orientation within the cavity, and includes terms 
that couple the rotational and translational motion. The 
vibrational-rotational Hamiltonian ""TiX^ is diagonal in 
the basis set \vJNLAMj\} with eigenvalues given by 
^£;y^ = hiuXiv + 1/2) + B^J{J + 1), where uj^ is the 
fundamental vibration frequency; B,j = Be — ae{v+l/2), 
where ae is an anharmonic correction to the rotational 
constant i?e[14]. Below 120 K, the thermally-activated 
rotational motion of the Cgo cages is suppressed jl3|, l28l | 
and ^V may be assumed to be time-independent. Ex- 
panded in multipoles it reads: 



"F(i?,rj) = 
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where the functions F are defined in Eq. [T] and ji takes 
even values. Terms with n = 2 constitute a harmonic 
potential energy function, while terms with n > 2 rep- 
resent anharmonic perturbations. Translation-rotation 
coupling terms are terms with non-zero l,j. 

All odd-j terms vanish for homonuclear diatomic 
molecules. For an icosahedral cavity, and assuming that 
longer-range intermolecular perturbations are negligible, 
all terms with odd A values vanish, as well as the terms 
with A = 2 and 4. We assume that all high-order terms 
starting from A = 6 are small and express the potential 
energy as ^V = ^V^ +^V' , where the isotropic harmonic 
term is given by "F° = "l^o'o^^oo + "^^o'o'^^^^oo and the 
perturbation due to translation-rotation and anharmonic 
coupling is given by "F' ^ ""V^o^R^F^g + ''Vo™4^4^oo 
The unperturbed Hamiltonian eigenvalues in the ba- 
sis \vJNLAMa) are given by E^j^^^j^,^^ 
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h^'u^iN + 1/2), where ^lo^ = {;"V^^^/i2Trm)) ' is the 
frequency for translational oscillations within the cavity. 
The matrix elements of "y' were evaluated analytically 
in the basis \vJNLAMa) using 100 states with N < 2 
and J ~ 1 , 3 for ort/io-H2 , and 60 states with N <2 and 
J = 0,2 for para-ii2- Matrix diagonalization leads to 
explicit but cumbersome expressions for the energy lev- 
els and eigenstates. A schematic energy level diagram is 
given in Fig[2l The ordering of the eigenvalues depends 
on the relative sign and magnitudes of the anharmonic 
term ^V^'o'^ and the translation-rotation coupling term 
"Kk?^- ^^^ ordering in Figl5]is consistent with the ex- 
perimental results. 
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FIG. 2: Diagram of H2 energy levels refined against low-T IR 
data. For the v — state only the ground state rotational and 
translational levels are shown. The energy levels of free H2 
are shown on the left while the effect of confinement by Ceo 
is shown on the right, both for para- and ortho-¥i.2- The ar- 
rows show transitions corresponding to the observed low-T IR 
peaks. The transitions marked by dashed lines are forbidden 
within the theory presented here. 



IR activity in H2@Ceo is due to a dipole moment, /i, in- 
duced by the constraining environment. The dipole mo- 
ment operator can be expanded in multipoles depending 
on the instantaneous H2 configuration [29[: 



t^q 
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where q denotes the spherical component and the A co- 
efficients describe the induced dipole moment. Since the 
dipole moment is a vector, there are restrictions on the 
allowed j and / values: (i) j + 1 must be odd, (ii) I = jztl 
from the triangle relation, (iii) for homonuclear molecules 
only even j terms are allowed. These restrictions imply 
selection rules for IR spectroscopy of H2@C6o- The selec- 
tion rule for the total angular momentum is AA = 0, ±1 
with the only allowed transitions having even values of 
A J and odd values of AL. In addition, because only the 
ground translational states (TV = L = 0) are populated 
at low temperature, the allowed transitions observed in 
the 6 K IR spectrum are to iV = L = 1 states (Fi g. ^ . 
The IR absorption amplitude at frequency w is |l4| : 

S^^Lu Y,P^inK,T)\{f\^lg\^)\^SiEf - E,~lu), (5) 
'■J 

where K ^ O or P selects ortho- or para-^2 ■ At the sam- 
ple temperature T, the fractional population pi(ni<-, T) of 
the initial state \i) is given by the Boltzmann distribution 
for ortho and para manifolds separately. Since the spin 
isomer interconversion is negligible for the endohedral 
complex [13|, the number of ortho and para molecules, 
no and np, is not in general governed by the Boltzmann 
distribution and must be determined empirically. Since 
all the observed transitions are from L = to L = 1, only 



TABLE L Experimental and calculated center frequencies, uj, 
and absorption line areas, St^ , of IR-active H2 modes at 6 K 
in H2@C6o. The forbidden transitions (Fig. la and dotted 
lines in Fig. 2) are used as frequency references for the fitting 
procedure. 
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terms with I — 1 and j — 0, 2 are to be considered in 
the dipole expansion. This implies that only A J = 0, +2 
transitions are observable. The dipole matrix elements 
for the observed transitions in H2@C6o between states 
|i) and |/) can be expressed as 



(/Im,IO 



An ^ 

^ J=0,2 



p'ixi^M, 



(6) 



where p^ = (*];'(r)^fi(i?)|AiJ(i?, r)|*^(r)*Jo(i?)) and 
[Xl^g]fi = {Fj:{Q)\Flf{n)\Fi{n)). The matrix elements 
reduce into a sum of products over radial integrals p' 
and known angular integrals. 

Three potential parameters Vo'o'^ ^^dc?^ ^Ko"^^ the 
rotational constant i?e, and p'^, p^^no/np were fitted to 
match the experimental frequencies and intensities. The 
potential parameters in the first vibrational state are 
{'■C': '■^00''} = {27 ± 6, 1.5 ± 0.2} J m-2 and ^V°°^ = 
(-2±20)102" J m-''. The low-T data is insufficient to de- 
rive accurately the anharmonic correction "'^Vpo'^ , which 
is poorly defined from the separation of A'^ = and 
iV = 1 levels only. This affects the value of ^Vj^g"^, 
because ^Vqq"^ and ^Vq^'^ are correlated. However, the 
translation-rotation coupling term ^Vqq^^ is well-defined. 
^V^Q^ and ^Vqo''* describe the potential within approx- 
imately 0.5 A from the Cgo cage center, which is the 
root square of the average square displacement for the 
A^ = 1 translational state. The fitted rotational constant 
is Be = 59.3 ± 0.2 cm"! while a^ = 2.98 ± 0.10 cm"! 
is obtained directly from the difference in the fundamen- 
tal vibrational frequencies for the ortho- and para-}l2. 
The ratio between the induced dipole moment param- 
eters is p^ / p^ — —2.0 ± 0.2. The ortho to para ra- 
tio no/np ~ 2.8 ± 0.2 is consistent with the equilibra- 
tion at any temperature warmer than 120K and suggests 
that there has been negligible spin isomer interconversion 
since the molecules were synthesized. The data and the 
best fit results are displayed in Fig. [1] and summarized in 
Table HI The frequencies of the pure vibrational transi- 
tions (Fig. [T^) are shifted by — 90cm~^ from the free II2 



value. The reduction in both the vibrational frequency 
and the rotational constant [30| are consistent with a 
predominantly attractive C-H interaction that slightly 
stretches the H-H bond. 

The results described here are in qualitative agree- 
ment with previous theoretical studies. The quantum- 
chemical calculation by Cross[8J gave a factor of two 
smaller translation-rotation coupling potential term and 
larger values for both harmonic and anharmonic poten- 
tial terms, leading to significant discrepancies in the IR 
line positions if used to reproduce the experimental data. 
Xu et aZ.|9[ used pair- wise Lennard- Jones potentials in 
their quantum-mechanical calculation. The ordering of 
the A sublevels is consistent with our results, even though 
the splittings are different. Xu et al. used two different 
potential energy surfaces, which gave different ^uiq , one 
smaller and the other larger than the experimentally ob- 
served ""luq ■ The splitting among the J = N = 1 energy 
levels, that is the measure of ^Vqq"^, remained larger than 
the experimentally observed value for both potential en- 
ergy surfaces. Although our measurements refer to the 
f — f excited state, while the numerical calculations are 
for the V = ground state, it seems that the theory us- 
ing a pairwise C-H potential is not accurate enough to 
describe the dynamics of H2@C6o. We believe that our 
results are a reference for theories modeling the interac- 
tion between H2 and curved carbon nano-surfaccs. 

Deformation of the cage, crystal field, or carbon iso- 
topomers in the cage may lower the symmetry and split 
the 4255 and 4261 cm~^ lines (Fig. [1)3) and cause the IR 
activity of the weak fundamental transitions (Fig. [T^) . 

In summary, the IR spectrum of endohedral H2 dis- 
plays a rich structure due to the coupled translational 
and rotational modes of the confined quantum rotor. 
Line positions and intensities are described by a theory 
involving multipole expansion of the confining potential 
and fitting of a small number of parameters. The next 
targets will be analyzing higher-T data to extract infor- 
mation about anharmonic corrections and the vibrational 
ground state, studying lower symmetry cages and differ- 
ent dihydrogen isotopomers. The accurate determination 
of the energy levels in the vibrational ground state may 
play a key role in understanding fully the low-T NMR- 
behavior of endohedral H2 in different fullerenes[13J. 
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